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Abstract We have analyzed "'^^^Cs decay data, ob- 
tained from a small sample onboard the MESSENGER 
spacecraft en route to Mercury, with the aim of setting 
limits on a possible correlation between nuclear decay 
rates and solar activity. Such a correlation has been 
suggested recently on the basis of data from ^'*Mn de- 
cay during the solar flare of 13 December 2006, and by 
indications of an annual and other periodic variations 
in the decay rates of "^^Si, "^^Cl, and ^^®Ra. Data from 
five measurements of the ^'^^Cs count rate over a period 
of approximately 5.4 years have been fit to a formula 
which accounts for the usual exponential decrease in 
count rate over time, along with the addition of a theo- 
retical solar contribution varying with MESSENGER- 
Sun separation. The indication of solar influence is then 
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characterized by a non-zero value of the calculated pa- 
rameter i, and we flnd i = (2.8 ±8.1) x lO'^ for i37Cs. 
A simulation of the increased data that can hypotheti- 
cally be expected following Mercury orbit insertion on 
18 March 2011 suggests that the anticipated improve- 
ment in the determination of ^ could reveal a non-zero 
value of ^ if present at a level consistent with other 
data. 
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In a recent series of papers (I Jenkins and Fischbach 



2OO9H Jenkins et al.lboOalFischbach et al.l2009HSturrock et al.^ 



2010at IJavorsek II et all bold ISturrock et al] bolObI ) 
evidence has been presented for a possible solar in- 
fluence on nuclear decay ra tes. Data analyzed by 
Jenkins an d Fisc hbachI (|2009l ) and lFischbach et al.l (|2009l) | 
indicate a possible correlation between the solar flare 
of 13 December 2006 and a decrease in the mea- 
sured decay rate of ^^Mn coincident in time with 
the flare. An analysis of data from Brookhaven 
National Laboratory (BNL) on the measured decay 
rates of ■^^Si and "^^Cl, and from the Physikalisch- 
Technische Bundesanstalt (PTB) in Germany on the 
measured decay rates of ^^^Ra and its daughters, 
show that both data sets exhibit similar annual vari- 



ations in their respec t ive de cay rates ([Jenkins et al 



2009t iFischbach etall boOOl). Similar periodic ef 



Falkenbertd (|2001l) . iBaurov et al.l (|2001 



fects ha ve been repo r ted by iParkhomov and Maklvaev 
(2 061: IParkhomovl (l2005l l2010b-G.d\ lElliJ dlQQOl) . 



2007), and 



citetshn98,shn00, and more recently by iJenkins et al 
( Submitted to NIM-A. 20 111) , in data from The Ohio 
State University. In addition to annual periodici- 
ties, evidence for other periodicities in decay data 
pos sibly associate d with so lar rotation is reporte d 



Sturrock et all (l2010al) . ISturrock et all (l2010bl) 



Fischbach et al.l (|201lD . and ISturrock et al.l (|201ll) . in- 
cluding evidence for a period of ^^33 days, and for a 2.11 



yr~^ Rieger-like periodicity. Since none of the rotation- 
related periodic signals (in what should be randomly 
distributed data) corresponds to any known terrestrial 
influence, these results support the inference of a so- 
lar origin to time- varying nuclear deca y rates, through 
some as yet unknown mechanism ([Fischbach et al 



2011 ) 



from which inferences can be drawn about the planet's 
surface composition. The Gamma-Ray Spectrometer 
(GRS) portion of the GRNS is a high-purity germanium 
(HPGe) sensor surrounded by a plastic scintillator anti- 
coincidence shield (ACS). During normal operation, the 
HPGe sensor is cooled to cryogenic temperatures using 
a miniature mechanical cooler. The mechanical cooler 
Th e interpretation of the data inlJenkins and Fischbachl has a limited operational life, with an expected mean 
( 20091 ). I Jenkins et al.l ( 20p9f). and 'Fischba ch et al. (200£ )| time to failure of approximately one year. Because the 

primary goal of the GRS is to map Mercury's compo- 
sition for one Earth year, the GRS has been used spar- 
ingly during MESSENGER'S seven-year cruise phase. 
The few times the GRS has been turned on have been 
to verify its operation and to make gamma-ray mea- 
surements of M ercury during the three Mercury flybys 
(lRhodesll201ll) . 



have b e en questioned b y Coopeii (|2009,) : iNorman et al 



(2008); ISemkow et al.l (|2009[ ). however, iJenkins et al 



20101 ) addressed all of those questions. Nonetheless, 



the nature of the effects reported in the original refer- 
ences remains uncertain at present, particularly since 
solar flares are not always corr elated in detectabl e 
changes in nuclear decay rates (|Parkhomovl l2010al) . 
However, if nuclear decays can in fact be influenced 
by the Sun, and specifically by the varying distance be- 
tween a decaying source and the Sun, then data from 
the MErcury Surface, Space ENvironment, GEochem- 
istry, and Ranging (MESS ENGER) mission to Mer- 
cury (jSolomon et al.l 120071 ) could lead to significant 
constraints on the magnitudes of such effects on the 
examined isotope, given that the MESSENGER-Sun 
distance has varied from 1.0689 AU to 0.30748 AU. The 
object of the present paper is to use decay data from a 
sample of "7Cs onboard MESSENGER to search for a 
possible variation of the ^'^^Cs deca y parameter A, or it s 



half-life Ti/2 = £n2/A = 30.07 yr ()Baum et al.l 120021 ). 



over the course of the mission to date. 

The MESSENGER spacecraft was launched on 3 Au- 
gust 2004, and Mercury orbit insertion (MOI) occurred 
as scheduled on 18 March 2011. The objective of MES- 
SENGER is to study Mercury from a highly eccentric 
orbit with a low-altitude periapsis. Following launch, 
several gravity-assist maneuvers were carried out by the 
MESSENGER spacecraft. These included a close flyby 
of the Earth approximately one year after launch (to 
reduce the required launch energy), and a subsequent 
deep-space maneuver to put MESSENGER on course 
for two Venus flybys to further lower the perihelion dis- 
tance. Three Mercury flybys (with a leveraging deep- 
space maneuver in between successive Mercury gravity 
assists) were used to gradually slow the spacecraft with 
respect to Mercury (A timehne for the MESSENGER 
mission trajectory is presented in Figure [ij) Among 
the seven instruments that comprise the MESSENGER 
physics payload is a Gamma- Ray and N eutron Spec- 
trometer (GRNS) (JGoldsten et al.l 120071 ). whose pur- 
pose is to map the elemental composition of the surface 
of Mercury. Galactic cosmic rays (GCRs) incident on 
the surface of Mercury produce characteristic gamma- 
rays and neutrons, which are detected by the GRNS, 
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Fig. 1 Timeline for the MESSENGER mission. Launch 
from Earth was on 4 March 2004, and Mercury orbit inser- 
tion (MOI) occurred on 18 March 2011. During the MES- 
SENGER mission the MESSENGER-Sun distance has var- 
ied from 1.0689 A.U. to 0.30748 A.U., a wide enough range 
to test the hypothesis of a time- varying nuclear decay rate 
that depends on distance from the Sun. Figure courtesy of 
NASA/ Johns Hopkins University Applied Physics Labora- 
tory/Carnegie Institution of Washington. 

Shortly after launch, readings from the GRS while 
in deep space revealed the presence of an unexpected 
small ^'^^Cs source. It is not known if this ^'^^Gs source 
is incorporated into the GRS sensor or some part of 
its housing. Nevertheless, because this source has been 
onboard MESSENGER from launch, its presence al- 
lows us, serendipitously, to set limits on the influence 
of solar activity on ^'^'^Cs decay. To date, flve sets of 
measurements of the ^'^^Cs decay rate have been made 
during periods denoted by Cruise through Cruise 4, 
and more are expected after MOI. In what follows, we 



develop a formalism to study a possible solar influence 
on ^■^^Cs decay, and we then apply this formalism to 
the existing data. We conclude with a discussion of the 
improvements that can be expected from data to be 
acquired after MOI. 



1 Theoretical Formalism 

The presence of the ^'^^Cs source onboard MESSEN- 
GER provides an opportunity to test whether the ^■^^Cs 
decay rate varies with MESSENGE R-Sun distance. 
Follow ing the discussion presented in iFischbach et al.l 
( 20091 ) we assume that the beta decay rate dN{t)/dt = 
N{t) of a sample containing N{t) unstable nuclei can 
be written in the form 



-dN{t) 
Jt 



A(t)7V(t) = [Ao + Ai(f,i)]A^(t), 



(1) 



where Aq represents the intrinsic contribution to the 
/3-decay rate arising from the conventional weak inter- 
action along with a possible time-independent back- 
ground arising from new interactions. The time- 
dependent perturbation Ai(r, i) is also a function of 
the distance r = jr'l between the decaying nucleus and 
the source, here assumed to be the Sun. For a pertur- 
bation whose influence on a decaying nucleus varies as 
1/r", where n is an integer, Ai(r, i) can be expressed 
in the form 



Ai(r,i) = AoC 



(n) 



R 



(2) 



where R = r(i — 0), and <^'^"-' is the parameter we are 
interested in constraining. In what follows we focus 
initially on ^'^^' = ^, which is the most likely variation 
that we expect, (i.e., an inverse-square law). Results 
for n = 1 and n = 3, which are the next most likely 
alternatives to n = 2, can be analyzed in a similar way. 
We note from Eqs. ([T]) and ([2]) that ^ is not a universal 
constant since it depends on the specified initial values 
R and i = 0. In the present paper we define t — as the 
launch time, in which case R = l AU = 1.495979 x 10^ 
km. Integrating Eq. ([Ij we find 



iV(i)=iVoexp{-Ao[t + e/W]}, 



(3) 



AoiVoexp(-Aot)<^l + e 



R' 

r2(t) 



dN{t) 
'df 



(4) 



We note from Eqs. p]|5)) that since /(O) = 0, the ef- 
fective ^■^^Cs decay constant at t == is A = Ao(l + ^), 
where A = 6.311(19) x 10~^ d^ is the standard value 
measured on Earth (at 1 AU) corresponding to T1/2 — 
30.07(9) yr. It is convenient to rewrite Eq. (j4]) in terms 
of directly measured quantities by first expressing Ao in 
terms of A, and then eliminating the unknown iVo by 
taking appropriate ratios. To lowest order in ^ <C 1 we 
find 



.At dN{t)/dt 
dN{t = Q)/dt 



= 1+e 



R' 



- XI{t) + Xt-l 



l + ^Bit). 



(6) 



We note from Eq. (O that since B(t = 0) = the 
count rate ratio approaches unity as i — > 0, as expected. 
Comparing Eqs. (|3]) and (jH]), the additional terms in 
the square brackets in Eq. ([6]) arise from computing 
the indicated ratio, and then expressing the final result 
in terms of the laboratory value A. 

In principle, the left hand side of Eq. ^ and B{t) 
are directly measurable, and hence can be used to set 
limits on ^. These limits can then be compared to 
those obtained in the laboratory. In practice, however, 
dN{t — 0)/dt was not measured, and hence our limits 
on ^ were determined from ratios of count rates ob- 
tained from measurements made during the five cruise 
periods. A summary of the data obtained from these 
measurements is presented in Table [TJ Among the en- 
tries in this table, the precision of the data which enter 
into the calculation of B{t) is such that the errors on 
B{t) can be neglected compared to those arising from 
the count rates (see discussion below). The integral 
I{t), which is exhibited in Figure [2l accounts for the 
cumulative ^-dependent contributions to the ^"^^Cs de- 
cay rate arising from both the change in MESSENGER- 
Sun distance, as well as from the normalization relation 

A = Ao(l + e). 

Returning to Eq. ([6]) , the experimentally interesting 
quantities are the ratios 

[e^p{XU)dN{U)/dt] 



[exp{Xto)dN{to)/dt]' 

where tg and ti [i = 1, 2, 3, 4) denote the start times of 
the five cruise measurement periods. These five mea- 
surements thus determine four ratios which, along with 
the corresponding values oi B{ti), can be used to obtain 



2 Analysis of the ^^^Cs Count Rates 



I{t) = / dt' 
Jo 



R 



r{t') 



(5) 



The MESSENGER GRS data analyzed in this paper 
are made publicly available in the NASA Planetary 



Table 1 Summary of Cs decay data from the MESSENGER mission. For each of the five measurement periods 
{to,ti,t2,t3,t4), these data include the calculated values of I{t), B{t), and A(t), for a perturbation varying as 1/r^ 



Cruise 
Initial checkout 



Cruise 1 
Mercury Flyby 1 



Cruise 2 
Mercury Flyby 2 



Cruise 3 
Mercury Flyby 3 



Cruise 4 
Final checicout 



Calendar date 


16 Nov 2004 


14 Jan 2008 


6 Oct 2008 


2 Oct 2009 


14 Apr 2010 


Julian date 


2453322.500 


2454478.500 


2454744.500 


2455102.500 


2455300.500 


Distance from the Sun (km) 


160599579.0277 


54595403.8012 


52751315.2631 


48132083.3548 


63523375.7902 


Elapsed time since Cruise (days) 


0.0 


1156.0 


1422.0 


1780.0 


1978.0 


Mean count rate in peak (cps) 


0.100±0.0012 


0.0942±0.0019 


0.0933±0.0020 


0.0886±0.0017 


0.0886±0.0020 


Total accumulation time (hours) 


60.00 


45.00 


40.00 


60.00 


39.00 


lit) (days) 


91.173 


2284.190 


3408.077 


5162.823 


6228.370 


B{t) 


-0.1316 


6.4397 


6.9235 


8.4530 


4.2842 


A(t) 


0.0000 


6.5751 


7.0551 


8.5847 


4.4159 
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Fig. 2 Plots of the integrals /(f), for n = 1,2,3, in Eq. 
|(2|. The integrals are defined in analogy to /"'(i) in Eq. 

Data System (PDS) archivqj according to the sched- 
ule given at the PDS web site. Table [1] lists the time 
periods for each of the five cruise measurements. Each 
of the data collection periods lasted for greater than 
24 hours, and some for over 60 hours. To ensure that 
the data analyzed here were acquired when the instru- 
ment parameters and background were reasonably con- 
stant, appropriate time windows were selected within 
each data collection period. For example, the Cruise 
period took place during the declining phase of a 
solar particle event, so the background was relatively 
high during the early portion of this collection period. 
Thus, the selected time window for the Cruise period 
was the last 60 hours of the collection period, when the 
background was the lowest. For the other data collec- 
tion periods, data were excluded during the Mercury 
flybys and when instrument parameters were varied for 
sensor testing. 

Figure [3] shows an example anticoincidence pulse- 
height spectrum from the Cruise data collection pe- 

^l|http://geo.pds.nasa.gov/missions/messenger/index.htm[l 



riod. Figure [5^ shows the electronic pulser peak used 
to correct for detector dead time; Figure [3}d shows the 
662 keV peak from the ^^^Cs decay. GRS dead time 
is largely driven by a 1000 Hz raw counting rate in 
the ACS. The dead time is monitored by a 7.62939 
Hz electronic pulse fed directly across the HPGe de- 
tector. The pulser line appears at a high energy in the 
pulse-height spectrum, far away from other prominent 
gamma-ray lines. The total counts in the pulser peak 
are obtained by fitting a linear function to the underly- 
ing background, subtracting this background from the 
pulse height spectra, and summing the net counts. The 
live time fraction is the ratio of the measured counts 
to the expected number of pulser counts. For the five 
measurement periods, the live time fraction varied from 
0.9 to 0.94. 

The total counts in the ^^^Cs peak are determined 
by dividing the spectra by the live time fraction, fitting 
a second-degree polynomial function to the background 
above and below the peak (dashed line in Figure [5}d) , 
subtracting this background, and then summing the net 
counts (dot-dashed line) within a window around the 
peak. The count rate uncertainties are calculated by 
assuming all uncertainties are due to Poisson counting 
statistics. Specifically, if the total net counts within 
the energy window are Njqet = Nspec - Nsack, where 
Nspec and Nsack are the total spectrum and back- 
ground counts within the energy window, respectively, 
then the total count uncertainty, a Net, is 



CTNet = \l(^Spec + ^ Back 



\/Nspec + NBack. (7) 

The corresponding count rate uncertainty is: 



(^dN/dt 



dN\ a Net 



dt I NNet 



(8) 



The uncertainty associated with the live time correction 
is neglected because its magnitude is almost an order 
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Fig. 3 Measured GRS pulse height spectra taken during 
the Cruise data collection period, (a) Pulser data (solid 
line) and background (dashed line) fit to the pulser data; 
(b) ^^^Cs data (solid line), background in the region around 
the ^^^Cs line (dashed line), and net ^"^^Cs spectrum (dot- 
dashed line). Solid vertical lines show the energy window 
over which the counts are determined. 



of magnitude lower than that of the uncertainty for the 
^^'^Cs peak. 

One other point regarding the GRS is that in or- 
der to reduce energy spreading effects from GCR ra- 
diation damage, the GRS is periodically annealed at 
high temperatures. Since launch, the GRS has been 
annealed four times (no annealing was carried out prior 
to Gruise 4). One consequence of annealing is that 
the active volume of the germanium crystal may shrink 
slightly, which in turn will slightly reduce the overall 
detector efficiency. By monitoring the position of the 
pulser peak after each annealing operation, we can es- 
timate the change in detector capacitance, and there- 
fore the active detector volume, by assuming a simple 
coaxial geometry for the HPGe detector. Since launch, 
the pulser peak has shifted ~7%, which equates to less 
than ~0.7% loss in detector volume; however, the ac- 
tual change in detector efficiency at 662 keV is not 



easily calculated (perhaps it could be better modeled 
through a Monte Garlo simulation) , but would likely lie 
between the estimates for the capacitance and volume. 
Such efficiency changes have not been included in this 
analysis. For future analyses, efficiency corrections can 
be monitored and/or included by observing changes in 
peak counts associated in time with annealing events, 
as well as by modelling the GRS efficiency response. 



3 Results 

Following the preceding discussion, we have determined 
£, by using Eq. ([6]) to write 



Mu^to) [dNiU)/dt] _ 
[dN{to)/dt] 



eA(f. 



1, 



(9) 



where i = 1, 2, 3, 4 and A{ti) = B{ti) - B{to) = B{ti) + 
0.1316. The starting dates U for Gruise 1, . . . Gruise 
4, are given in Table [1] as are the calculated values for 
the integrals I{ti) in Eq. (P, and A{U) in Eq. (0. 
Gombining these results with the four ratios obtained 
from the measured countrates allows us to fit Eq. ([9]) 
to a straight line, yi = £,Xi + b, whose slope is ^. We 
find, 

^ = C^2^ = (2.8 ± 8.1) X 10"^; 6^2) = 1.02 ± 0.03 (10) 

The result for ^ in Eq. (|10p can be compared to 
the annual variation in the '^^Si/^^Cl coun t rate ratio 
reported by iJenkins arid Fischbach ! f2009) from their 
analy sis of the data of [ Alburger et al. (1986 ) taken at 
BNL. Uenkins and FischbachI ( 2009 ) found a fractional 
change of ~ 3 x lO'^ in N{^'^Si)/N{^^Gr) from peri- 
helion to aphelion, and this level of sensitivity may be 
achievable following MOI given the result in Eq. (fTO|) 
as we now discuss. 



4 Post-MOI Simulation 



The result in Eq. (|10|) suggests that the availability 
of a much larger data set following MOI on 18 March 
2011 could lead to significantly improved constraints on 
£_ for ^'^^Gs. In this section, we simulate post-MOI data 
that can be expected in order to estimate the level of 
improvement on ^ that might be realized. 

We can estimate the post-MOI countrates by re- 
expressing Eq. ([6]) and (HJ in the form 



XAt 



dN{ti)/dt 



1 + C 



i?2 



r^iU) 



A / dt' 



R 



it') 



_dN(ti)/dt 
\At-l-D 



(11) 



where 



D = 



R^ 



'{h 



-X 



dt' 



R 



r{t') 



Xti -1 = 9.1542,(12) 



accounts for the accumulated change in the ^'^^Cs 
count rate between Earth launch at i = and MOI at 
t = ti. In Eq. (fTTj) . ti is the time of the i*'' measure- 
ment, At = ti — ti in days since MOI, and we have used 
for r (t) the actual (planned) trajectory ephemeris, cal- 
culated by the Jet Propulsion Laboratory HORIZONS 
system. The simulation was run for slightly more than 
one Earth-year following MOI, during which each mea- 
surement was assumed to last for 11 days. We have 
also assumed th at £ = —0.0031. whic h is compatible 



with the data of lAlburger et al.l (|l986f ). and have fixed 



the initial ^'^^Cs count rate to be 0.08 cps, which is sug- 
gested by the data in Table[T] Under these assumptions, 
the standard fractional statistical error in each 11-day 
run is approximately ANi/Ni = 0.0036, and hence the 
fractional error in the ratio given in Eq. (|lip is approx- 
imately V2(0.0036) = 0.0051. A uniformly distributed 
random error was then added to each expected mea- 
surement to simulate uncertainties in the observations. 
The objective of the simulation is to recover the input 
value S, = —0.0031, and to then determine its associated 
uncertainty. 

The results of our simulation are shown in Figs. S] 
andEl In Fig. HI the error bars represent the 1/Vn 
errors in each point, and the asterisks (*) denote the 
result of a random measurement of the count rate ratio 
within the I/^/N range. The solid line represents the 
expected results in the null-case (^ = 0), correspond- 
ing to a pure exponential decay in time. Fitting the 
simulated data to Eq. (jTTj) . we find (see Fig. [5]), 



C= (-3.1 ±0.4) X 10" 



(13) 



for the post-MOI data that we have generated. We 
note from Eq. (fT3|) that the central value of ^ is in 
agreement with the input value ^ — —0.0031 that we 
have used. More significantly, the uncertainty in ^ is 
sufficiently small to suggest that a non-zero value of ^ 
near the nominal level that we have assumed might be 
detectable using a more sophisticated analysis of the 
actual post-MOI MESSENGER data. 
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Fig. 4 Simulated count rates as a function of time. The 
error bars represent the 1/vN statistical uncertainties at 
each point, and the asterisks (*) are the results of simulated 
"observations". The solid line is the prediction for ^ = 0, 
corresponding to a pure exponential decay in time. 
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Fig. 5 Determination of ^ from the simulated data. The 
vertical and horizontal axes correspond, respectively, to the 
left-hand and right-hand sides of Eq. (jlip . where d = Ni, 
etc., and the solid line gives the best fit to the simulated 
data. The inferred value of ^ is given in Eq. (jlSp . 



References 

Alburger, D.E., Harbottle, G., Norton, E.F.: Earth and 
Planetary Science Letters 78(2-3), 168 (1986) 

Baum, E.M., Knox, H.D., Miller, T.R.: Nuclides and Iso- 
topes: Chart of the Nuclides, 16th edn. Lockheed Martin 
(KAPL), (2002) 

Baurov, Y.A., Konradov, A. A., Kushniruk, V.F., 

Kuznetsov, E.A., Sobolev, Y.G., Ryabov, Y.V., Senke- 
vich, A. P., Zadorozsny, S.V.: Modern Physics Letters A 
16(32), 2089 (2001) 

Baurov, Y.A., Sobolev, Y.G., Ryabov, Y.V., Kushniruk, 
V.F.: Physics of Atomic Nuclei 70(11), 1825 (2007) 

Cooper, P.S.: Astroparticle Physics 31(4), 267 (2009) 

Ellis, K.J.: Physics in Medicine and Biology 35(8), 1079 
(1990) 

Falkenberg, E.D.: Apeiron 8(2), 14 (2001) 

Fischbach, E., Jenkins, J.H., Sturrock, P.A.: arXiv 
(arXiv:1106.1470vl [nucl-ex]), 1 (2011) 

Fischbach, E., Buncher, J., Gruenwald, J., Jenkins, J., 
Krause, D., Mattes, J., Newport, J.: Space Science Re- 
views 145(3), 285 (2009) 

Fischbach, E., Sturrock, P.A., Jenkins, J.H., Javorsek II, 
D., Buncher, J.B., Gruenwald, J.T.: In: Kostelecky, 
V.A. (ed.) Proceedings of the Fifth Meeting on CPT and 
Lorentz Symmetry, p. 168. World Scientific: Singapore, 
(2011) 

Goldsten, J., Rhodes, E., Boynton, W., Feldman, W., 
Lawrence, D., Trombka, J., Smith, D., Evans, L., White, 
J., Madden, N., Berg, P., Murphy, G., Gurnee, R., Stro- 
hbehn, K., Williams, B., Schaefer, E., Monaco, C, Cork, 
C, Del Eckels, J., Miller, W., Burks, M., Hagler, L., De- 
Teresa, S., Witte, M.: Space Science Reviews 131(1), 339 
(2007) 

Javorsek II, D., Sturrock, P.A., Lasenby, R.N., Lasenby, 
A.N., Buncher, J.B., Fischbach, E., Gruenwald, J.T., 
Hoft, A.W., Horan, T.J., Jenkins, J.H., Kerford, J.L., 
Lee, R.H., Longman, A., Mattes, J.J., Morreale, B.L., 
Morris, D.B., Mudry, R.N., Newport, J.R., O'Keefe, D., 
Petrelh, M.A., Silver, M.A., Stewart, C.A., Terry, B.: As- 
troparticle Physics 34(3), 173 (2010) 

Jenkins, J.H., Fischbach, E.: Astroparticle Physics 31(6), 
407 (2009) 

Jenkins, J.H., Mundy, D.W., Fischbach, E.: Nuclear Instru- 
ments and Methods in Physics Research Section A: Accel- 
erators, Spectrometers, Detectors and Associated Equip- 
ment 620(2-3), 332 (2010) 

Jenkins, J.H., Fischbach, E., Buncher, J.B., Gruenwald, 
J.T., Krause, D.E., Mattes, J.J.: Astroparticle Physics 
32(1), 42 (2009) 

Jenkins, J.H., Herminghuysen, K.R., Blue, T.E., Fischbach, 
E., Javorsek II, D., Kauffman, A.C., Mundy, D.W., Stur- 
rock, P.A., Talnagi, J. A.: Additional experimental evi- 
dence for a solar influence on nuclear decay rates (Sub- 
mitted to NIM-A, 2011) 

Norman, E.B., Browne, E., Shugart, H.A., Joshi, T.H., Fire- 
stone, R.B.: Astroparticle Physics 31(2), 135 (2009) 

Parkhomov, A.G.: International Journal of Pure and Ap- 
plied Physics 1(2), 119 ( 2005) 

Parkhomov, A.G.: ArXiv (|arXiv: 1006.22951 /1 [physics.gen- 
ph]), 1 (2010a) 



Parkhomov, A.G.: ArXiv (ar Xiv: 1010. 1591 vl [physics. gen- 

ph]) (2010b) 
Parkhomov, A.G.: ArXiv (.arXiv:1012.4174 vl [physics.gen- 

ph]) (2010c) 

Parkhomov, A.G.: ArXiv (ar Xiv:1004.1761 vl [physics.gen- 

ph]) (2010d) 
Parkhomov, A.G., Maklyaev, E.F.: Fizicheskaya Mysl 

Rossu (1), 1 (2004) 
Rhodes, E.A.: In preparation (2011) 
Semkow, T.M., Haines, D.K., Beach, S.E., Kilpatrick, B.J., 

Khan, A.J., O'Brien, K.: Physics Letters B 675(5), 415 

(2009) 
Solomon, S., McNutt, R., Gold, R., Domingue, D.: Space 

Science Reviews 131(1), 3 (2007) 
Sturrock, P.A., Buncher, J.B., Fischbach, E., Gruenwald, 

J.T., Javorsek II, D., Jenkins, J.H., Lee, R.H., Mattes, 

J.J., Newport, J.R.: Astroparticle Physics 34(2), 121 

(2010a) 
Sturrock, P.A., Fischbach, E., Jenkins, J.H.: Solar Physics 

DOI: 10.1007/sll207-011-9807-5 (2011) 
Sturrock, P.A., Buncher, J.B., Fischbach, E., Gruenwald, 

J.T., Javorsek II, D., Jenkins, J.H., Lee, R.H., Mattes, 

J.J., Newport, J.R.: Solar Physics 267(2), 251 (2010b) 



This manuscript was prepared with the AAS lATgX macros v5.2. 



